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Abstract. We shall discuss cosmological models in extended theories of gravitation. We 
shall define a surface, called the model surface, in the space of observable parameters which 
characterises families of theories. We also show how this surface can be used to compare 
with observations. The model surface can potentially be used to falsify whole families of 
models instead reasoning on a single model basis as it is usually done by best fit arguments 
with observations. 
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1 Introduction 

The aim of this paper is to provide a detailed framework for interpreting extended cosmologies 
in the family of Palatini /( 7 ^)-models. The interpretation of extended theories is trickier than 
usually recognised in the literature, since the dynamics is not sufficient to constrain all details of 
the interpretation or the relation between the cosmographic parameters and the corresponding 
model quantities. All interpretation issues are generated by the fact that in Palatini f{IZ)- 
theories one has two metrics which play a natural role, the original metric g and a conformal 
metric g which naturally appears when solving the model. 

Prom the geometric point of view, while in standard general relativity (GR) the geometry 
of spacetime is described by a Lorentzian structure, in /( 7 ^)-theories spacetime is described by 
a Weyl geometry, i.e., by a conformal structure [s'] together with a compatible connection E 
which, in view of field equations, turns out to be the Levi-Civita connection of the conformal 

*This paper is a contribution to the Special Issue on Analytical Mechanics and Differential Geometry in honour 
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metric g= ip ■ g. The scalar field ip is called the conformal factor and it parametrises conformal 
transformations. Weyl geometries have been considered in the past in gravitational theories; 
see, for example, [4, 23, 24, 26, 27, 28]. 

Unlike in standard GR, where one has only one metric which does everything, in Palatini 
/(77)-theories one has to declare which metric is used to do what; see, for example, [1, 2, 7]. 
For example, the dynamics of field equations does not determine the free fall. If one couples to 
a matter field of course one can compute the characteristics of matter field equations and argue 
that characteristic curves describe free fall. In particular, if one couples a Klein-Gordon scalar 
field (p to the metric g then characteristics are determined as the geodesics of g. Unlike lightlike 
geodesics, the timelike geodesics of g are not geodesics of g. Accordingly, it seems that the free 
fall is determined by the variational principle assumed for gravity and matter. 

However, one can show that a conformal transformation can be extended to act on the 
matter field cp as well and define a new matter field cp = tp'^ ■ cp (with the parameter a G M to 
be determined) so that the new field (p obeys new (different) field equations (since conformal 
transformations are not symmetries of the matter field equations). The exponent a in the 
action of the conformal transformation on the Klein-Gordon field can be chosen so that the 
characteristics of the new field cp are geodesics of the conformal metric g] see [18]. Something 
similar can be done for Dirac fields. 

Hence arguing that the variational principle determines the free fall is correct, provided one 
assumes what is observed as a matter field. Depending on the choice oi <p oi <p as the observed 
field, then free fall is determined to be described by g or g, respectively. In this paper we shall 
assume that free fall is given by g. 

Other examples are the issues of whether the dynamics determines the metric which is used for 
defining spatial lengths, minimal coupling to matter, standard clocks. For example, in standard 
GR it is natural assuming that physical lengths are associated to the geometric length measured 
with respect to the only metric one has around, namely g. Having more that one metric in the 
model makes the choice much more problematic. 

Let us consider two parallel mirrors at a constant distance and a light ray bouncing between 
them, which is called a gravitational clock. If one considers an atomic clock (which is in fact 
built with some quantum harmonic oscillator) and compares its rate to the rate of a gravitational 
clock, the assumption of standard GR is pretty neat: the two clocks (at the same point) are 
synchronised forever. 

It is quite obvious that the discrepancy between the two clocks, if any, cannot be too big other¬ 
wise we shall see some effect at the solar system scales. However, there is no strong constraint 
on what happens at bigger scales. One should not assume that there is no discrepancy. On the 
contrary, the correct procedure would be to assume there is, derive observational consequences, 
and constrain the effects by observations. 

We shall hereafter follow this second guideline. In this paper we shall assume that lengths are 
related to g, not to g. Accordingly, we shall assume that matter is coupled to g, and not to F. 

These choices are also motivated by the Ehlers-Pirani-Schild (EPS) framework for gravita¬ 
tional theories; see [14]. In the early seventies EPS proposed a framework able to derive the 
geometry of spacetime based on basic properties of light rays and (free falling) particles. In this 
framework one does not assume a Lorentzian structure on spacetime, rather a Lorentzian metric 
(actually a Lorentzian conformal structure) on spacetime can be defined due to the properties 
assumed for light rays. As a consequence, in the EPS framework there is a tight link between 
the geometric structure of spacetime and observations. One could say that the EPS frame¬ 
work instructs an observer to establish observational protocols and to define distances and time 
lapses. Once this is done EPS considered the properties of particles and find they are described 
by a connection F (actually by a projective structure on the spacetime) which describes, by 
construction, the free fall. 
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The two structures ([ff], f) on spacetime are not arbitrary though. Since one knows that light¬ 
like directions are an upper bound for the speed of particles, one can show that the connection T 
must be in the form 

n, = {gVp, + \{9^^9P, - ( 1 - 1 ) 

for some covector A^. This is quite similar to the original proposal by Weyl for a unified theory 
of gravity and electromagnetism (then proven wrong); see [28]. However, we are not interpreting 
here A as an electromagnetic potential and for us the covector will be related to the possible 
shift in synchronisation between atomic and gravitational clocks. 

On the basis of the EPS framework we shall consider gravitational theories for the fields (^, T) 
to be extended theories of gravitation (ETC) when their field equations enforce as a consequence 
the relation (1.1); see [10, 11]. We shall call such theories integrable when field equations also 
imply that A is an exact 1-form and as a consequence T is a metric connection; see [13, 16]. 

In principle one can consider any invariant of the curvature to define a dynamics. We shall 
hereafter restrict to Palatini /(7^)-theories in which the dynamics is given in terms of a (regular) 
function of the scalar curvature TZ. One can show that Palatini /(7^)-gravity is an integrable 
extended theory of gravitation (iETG). In view of the EPS framework we shall interpret T as 
defining the free fall, while g is related to distances and time lapses. 

Standard GR is then a particular case of Palatini /(7^)-theory. However, the standard 
Hilbert-Einstein dynamics imposes that the conformal factor ip is exactly and everywhere equal 
to 1. In this special case the two metrics g and g do coincide and all issues related to choosing 
which metric is used to do what are overcome. 

It is particularly worth noticing that the EPS project was originally started to show that 
standard GR follows from assumptions and it almost reached the goal. Instead of finding 
a Lorentzian structure on the spacetime as assumed in standard GR, EPS find something less 
(just a conformal structure of Lorentzian metrics) and something more (a projective structure 
of connections which is not uniquely determined as the Levi-Givita connection of the metric 
structure). 

EPS asked in their analysis what one needs to add to the framework in order to force the 
structure to be Lorentzian. They showed that in fact if one requires that atomic and gravitational 
clocks have the same rate then standard GR follows. 

In view of this, Palatini /(7^)-theories can be seen to generalise standard GR in two respects. 
In the first place, the modified dynamics induces (effective) dark sources which can be used 
to model dark matter and dark energy. In the second place, the conformal factor does account 
for the possibility that gravitational clocks shift with respect to atomic clocks, the effect being 
measured by the field (p, i.e., the conformal factor; see [9, 17]. 

These two effects (dark sources and time shifts) are quite independently assumed in the 
theory and in our opinion one should not assume they are not there. Instead, they should be 
analysed by considering them in a model, one should make predictions and test the results, 
possibly ruling them out. 

The material is organised as follows. In Section 2 we shall review Palatini /(77)-theories. The 
conservation of stress tensors and their relation with conformal transformations is discussed in 
Appendix A. In Section 3 we shall discuss cosmological models based on Palatini /(77)-theories. 
In Section 4 we shall consider the example f{TZ) =TZ — §77^, the solutions of which are briefly 
discussed in the Appendix B. 

2 Palatini /(7?.)-theories of gravitation 

In view of the EPS formalism one should consider a Palatini formulation. For the sake of 
generality, one can consider a general dynamics, not restricting a priori to standard GR. For 
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the sake of simplicity let us consider a Palatini /(7^)-dynamics, and we shall assume for it the 
interpretation described above. 

Thus our theory is described by two fields (fi'/xi/, F^^) a metric g which is associated with rulers 
and clocks and a (torsionless) connection F which is associated to free fall of test particles. 

Let us define the scalar curvature TZ = and the Lagrangian 

L = f{n)y/gda + L^at, 

where da denotes the standard local pointwise basis of 4-forms on a 4-dimensional spacetime M 
induced by coordinates. The field equations turn out to be 

- y{T^)9iMu = T^u, 

= 0 , 

where the stress tensor is defined by variation of matter Lagrangian with respect to the 
metric g. We also assumed that the matter does not couple to the connection F but only to the 
metric g and we understand the matter field equations. 

The second equation does determine the connection F as the Levi-Civita connection of a con¬ 
formal metric 

9fiu ■ — T — l^l 

for some constant (/^o chosen so that one has ^{xq) = 1 here and now. Let us assume that the 
conformal factor (p is positive and the metric g is also Lorentzian. Using this information in 
the first equation, the Ricci tensor becomes the (symmetric) Ricci tensor of the conformal 
metric g, i.e., one obtains 

nn)R^y -\f{n)g^y = V. 

In order to solve this equation let us consider its trace by g^'^ first 

/'(7^)7^-2/(7^) = ^, (2.1) 

which is called the master equation. It is an algebraic (i.e., not differential) equation in TZ and 
T = gi^'^T^y-, except some degenerate cases, for a generic (analytic) f(TZ) one can solve it for TZ 
obtaining 

TZ = TZ{T). 

Accordingly, one can express the conformal factor as a function of the trace T, i.e., 

f'injT)) 

Using again this information to write the first field equation in terms of the new conformal 
metric g one has 

Rfiu ~ T 2^f{'TZ)g^y — i^TZgfj^y, 

where R = (p~^TZ is the scalar curvature of the conformal metric. In other words, the conformal 
metric g is determined by a Einstein equation though with a modified source stress tensor 

•= T ~ (2-2) 

which can in fact be considered as a function of g and the matter fields (included the conformal 
factor p). 

As a minimal check one can show that if the matter energy-momentum stress tensor T^y is 
conserved (with respect to g) then the effective stress tensor T^y is automatically conserved as 
well, with respect to g, though. See Appendix A for details. 
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3 Palatini /(7?-)-cosmologies 


We can use this theory of gravitation for producing cosmological models. For, one has to state 
cosmological principle, i.e., restricting to spacetimes with an isometry group which enforces 
spatial homogeneity and isotropy (or, equivalently, spatial isotropy with respect to any point). 
As Benenti pointed out, the cosmological principle freezes dynamics since there is only one 
dynamics, namely the one dictated by Friedman-Robertson-Walker (FRW) equations, which is 
compatible with homogeneity and isotropy; see [5]. 

However, if, as in this case, a theory can be described in terms of more than one metric, one 
has to state for which metric isometries are considered in the cosmological principle. Luckily, in 
Palatini /(7?,)-theories this is not an issue, since the conformal factor (p(t) itself is expected to be 
a function of time only, there are two metrics which are conformal to each other by a conformal 
factor which depends on time only, and if a metric g is spatially homogeneous and isotropic then 
the conformal metric g is spatially homogeneous and isotropic as well, though in a particular 
and precise sense which we need to make explicit for later convenience. 

A metric g obeys cosmological principle if and only if there exists a coordinate system 
{t,r,6,4>) in which the metric is in the form 

+ 


dVi^ := dO"^ + sin^(0)d(/>^ 


(3.1) 


g = —dt^ + a^{t) 


dr^ 


1 — kr^ 


for some function of time a{t) which is called the scale factor. A metric in this form is called 
a FRW metric. One can (and we will) always rescale the radial coordinate so to have a(to) = 1 
at a time to denoting today (which for convenience can be set to to = 0). 

If the conformal factor is a function of time only the conformal transformation sends a FRW 
metric into a FRW metric, though with another time and scale factor. If g is in the form (3.1) 
with radial coordinate rescaled to obtain a(to) = 1, then the conformal metric g = (p{t) • s' is in 
the form 


g = —dP + cP{t) 


dr'^ 

1 — kr"^ 



ait) = vV(t)a(t), 
di = ip{t)dt, 


(3.2) 


and one still has d(to) = 1 if and only if the conformal factor (^(t) is chosen so that ipito) = 1. 

Most of what follows relies on having a metric in FRW form and, consequently, it applies 
to both g and g. We shall hereafter briefly review what is routinely done for g in standard 
cosmology with the aim of then applying it to g. 

Once one has a FRW metric one can compute its Einstein tensor and find it in the 
form of a stress tensor of a perfect fluid i.e., there exists a vector field u which is a timelike 
unit vector for the metric g orthogonal to the submanifolds at t = const and two functions p(t) 
(the pressure) and pit) (the energy density) such that 


Ti,u = ip + p)u^Uy + pg^y 


T = -{p + p) + Ap = 2,p - p. 


(3.3) 


This has not much to do with dynamics since the stress tensors for a perfect fluid are the 
only tensors which are compatible with spatial homogeneity and isotropy. In some sense the 
cosmological principle freezes the dynamics and reduces it to a purely kinematical fact: any 
FRW metric obeys Einstein equations with a perfect fluid source for a suitable pressure and 
energy density functions. The dynamics only specifies the relation of density and pressure with 
the matter fields. 

Einstein equations are 10 equations, though in this case only two of them are independent. 
They are the so-called Friedman equations 


d\ kc^ 1 
o / 3 


d 1 , 

- = -x(/3 + 3p . 
a 0 


(3.4) 
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These are two equations for the three functions a{t), p{t), p{t) which then cannot be deter¬ 
mined unless a further equation is provided. In principle one should obtain that extra equation 
from the matter field equations. In practice, it is customary in cosmology to describe matter 
by means of an (algebraic) equation of state (EoS) which constrains p and p. The most general 
EoS a{p, p) = 0 usually is chosen in the form p = p{p) and often linearised simply to p = wp, 
where w is called the parameter of the EoS, the value of which determines the kind of matter 
we are dealing with. Different cases will be analysed below. 

Equations (3.4) are Lagrange equation for a suitable Lagrangian. If we consider the mat¬ 
ter described by an energy-density function p{a), then we can define the so-called point-like 
Lagrangian 


L = (—od^ -|- kc^a — Xp{a)a^')dt, 


(3.5) 


which is obtained by integrating on space the Hilbert Lagrangian for a ERW metric (and neglec¬ 
ting a total time-derivative). 

This Lagrangian is time-independent so that its “total energy” 


Ti = —2aa^ -|- aa^ — kc^a -|- Xp{a)a^ 


of kc^ , / N t 
^ H- 5 — Xp[a) o' 


is conserved. Notice that by imposing this quantity to be zero one obtains 

cf kc^ 

“2 -T “ 

which agrees with the first FRW equation by setting A = g. The equation of motion for the 
Lagrangian (3.5) is 

2d = l{2p + pa)a => - = ^{2p + pa), 

3 a D 

which becomes the second FRW equation provided one has the pressure defined to be 

-{p-\-3p) = 2p-\- pa p = - [p-\- ^pa). 

This condition is equivalent to the conservation of the energy-momentum stress tensor, namely 
= 0 for a fluid tensor in the form (3.3). 

Thus we select the type of matter fluid by choosing the function p(a) which in turn determines 
the matter equation of state (EoS) p = p{p). The definition of the pressure above is in fact 
a form for the EoS parametrised by the scale factor a 


p = p{a), p = p{a) ■.=-{p{a) + \p{a)a). 

Eor example, dust is determined by setting p = 0, i.e., 3p p'a = 0 from which p = p^ai 
follows. More generally, assuming an EoS in the form p = wp one has 


P — Poa 


3(io-l-l) 


Radiation is defined to be p = ^p, thence p = poa~^. Cosmological constant or {quintessence) is 
defined to be p = —p, so that one has p = po; see [9]. 

One can also assume the so-called polytropic equation of state p = p^ and obtain 


p{a) = 


3 


1 

7-1 


(ca)3(7-i) - 3 
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In view of the point-like Lagrangian, FRW cosmology is captured by first FRW equation 
which is a Weierstrass equation: 

^ ^ — \p{a) = 0 = —k + \p{a)a^ =: d>(a). (3-6) 

This function d>(a) is essentially the so-called effective potential. The quantity 



a 


is called the Hubble parameter. The value of the Hubble parameter today is denoted by Hq. 

One can observe that if A: = 0 in equation (3.6) then p'f = 377q which is called the critical 
density (today). Hence if one has different matter species each with its own density pi = p^oT'' 
and pressure pi, and its own EoS pi = Pi{pi), the equation (3.6) specifies to 


= o\ P 


3k 



Pi + PO j — gP, 


where we introduced one more species, called the curvature species, with density pk = —3ka~'^ 

k 

The total density p = Pi ™ust be critical today. 
i=0 

Then one can recast the FRW equation as 


5 = 5P5'E5 = 'ff»E«?“"-. 

Po 


i=0 


where we set 


= — 
cr ’ 

Pa 




i=0 


which are called abundances (today) and are one-dimensional parameter for each species. Since 
a set of abundances are constrained by X] they can be interpreted as percentages of 

each species in the total budget corresponding to the critical density. A set of such percentages 
are called a cosmic pie. Then in this more general framework, the FRW equation reads as 


= a^Ho = $(a; Hq, 

i 

and solutions are uniquely determined by initial condition a( 0 ) = 1 and today abundances H?. 


3.1 Einstein frame 

In Palatini /(77.)-cosmologies, also the conformal metric 5 is a FRW metric and as such it obeys 
Friedman equations for some suitable pressure and density as well. We already showed that the 
conformal metric can be recast in FRW form by a suitable rescaling. 

Thus we have the metric in the from (3.2). Since it is a FRW metric, it obeys some sort of 
Friedman equations which means that the modified stress tensor must be in the form of 
a perfect fluid stress tensor. 

In fact the timelike unit vector u for the metric g is 


u^g^uu'' = -1 




u^gfjiuu'' = -1. 
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The corresponding covector is 


Wm = dfiuU"' = ipg^j.uu" = = ^/^u^ 


and the modified stress tensor (2.2) is in the form 


+ -^{p + p - ifTVjg^iv, 


ip- 


Aip^ 


where the master equation has been used to eliminate the dependence on f{TZ). 

Thus also is in canonical form provided one defines the effective pressure p and energy 
density p as 


P + P 


P + P 
(^2 


3(p + /o) ~_P + P 

4 (/?’ ^ 4 (/?^ Aip 


By replacing TZ = TZ{T) = TZ{T) = Tl{3p — p) and p = ip^ ^f'{TZ{T)) = ^f'{TZ{3p — p)) and 

using the EoS, one obtains the effective pressure and energy density as functions of p, i.e., 


p = Kp), p = p{p), 

which is an implicit form of the EoS for the effective matter (parametrised by the matter energy 
density p). 


4 An example 

Let us consider the simple case (see [22]) 
f{R) = n- e>0. 

The stress tensor is assumed to be a fluid stress tensor. Visible matter will be described 
by this tensor and, for the sake of simplicity, we will assume visible matter to be dust. The 
master equation reads as 

7Z — eTZ^ — 2TZ + eTZ^ = —TZ = —p =► TZ = p 


and the conformal factor is 

1 — eTZ 1 — ep 

P = -= , 

(po 1 - epo 


P — doo 


Let us stress that there are two allowed regions for the scale factor, one for a > am = (epo)^^^ > 0 
and one for negative values of a. The effective pressure and energy density are 


P = 

P = 


3 + (^ 

Aip'^ 

1 - P 

4(y92 


1 - epo 
(l-ep)2 

e 1 - epo 
4 (1 - ep) 


(1 - fepo - 

^{p- po)p. 


hp)p, 


This corresponds to an effective EoS; see Eig. la. 

The relation between the scale factors is (see Eig. lb) 


d(a) 




1 - ep 
1 - epo 


a. 















Extended Cosmologies 


9 




Figure 1. a) Effective EoS p{p). (The dash line corresponds to a negative energy density p for visible 
matter, which may be considered unphysical.) b) Scale factors a(a). Dash line is standard GR. 


At this point we can expand p and p as functions of a and check that 


P=-[P+\%~^ 


dp dp dp da 
dd dp da da 


-3poa 


-4 


dp 

dp 



(4.1) 


The condition (4.1) is in fact generally equivalent to conservation laws for the effective energy- 
momentum stress tensor which is general follows from conservation of the original energy- 
momentum stress tensor T^y (see Appendix A), which in turn follows from EoS imposed for 
matter, or, equivalently, from covariance of the matter Lagrangian. 

Then one decides what kind of effective matter is needed in the model. This is not something 
decided within cosmology, it can be a choice of ours or an input from other physics (e.g., 
constraints from baryogenesis or structures formation). 

Eor example, let us assume we want a model with cosmological constant, dust and radiation. 
This corresponds to a total density in the form 


P'd 


Pr 




with the parameters {p^, p^, pP) to be determined so that p best approximates p. 

Best fit can be obtained simply by imposing that the value, first and second derivatives of p 
and p coincides at a = 1, i.e., today. The best approximation takes the form of 


Pa = Pa{(^:Po), = pJ](e,po), p° = p°(e,po)- 


(4.2) 


This functions are quite complicated functions, also for a simple model as the one we are 
considering here. However, one can have an idea of their shape (and the effective physics they 
describe) by expanding them for small e (i.e., small modification of standard GR). In this case 
we have: 


pi- 



^2 
4 ’ 




Pr 




9Po _ 9Po 2 
8 8 


Let us remark that those reduce to visible matter (dust only) for e = 0 and that all energy 
densities are positive for a small positive e. 

This means that Palatini /(7^)-cosmology with f{Tl) = TZ — and purely dust visible 
matter is well approximated today by standard cosmology with cosmological constant, dust and 
radiation. Of course, far in the past and future the approximation fails to be accurate. 
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Figure 2. Line model for f{Tl) = TZ — §7^^- Dash line is the constraint — 1- The dot is the 

observered cosmic pie (e.g., WMAP). 


We can also compute the abundances of different kinds of effective matter by dividing by the 
critical density today p" = 3Hq = po to obtain the effective abundance parameters today 


^“8 4 " ’ 




1 - 


5po 


e + 


IlPo 2 

^ t 1 




9po 


-e — 


9 Po 2 
8 


In principle, one can try and determine e so that, e.g., Ha — 0.75 as observed, e.g., by WMAP. 
However, there are numerical evidences that the (exact) abundance of effective dark energy in 
f{TZ) = TZ — models cannot exceed a value of about 0.62. In this case the whole family of 
models f{TZ) =Tl — can be rejected since it is not compatible with observations. 


4.1 The model surface 

What we got above as a description of the cosmological model is the abundance parameters 
(HO, HO, HO, HO,...) as functions of the parameter e which in fact parametrises the family of 
models under consideration, namely in the simple example of the function f(TZ) = TZ — §77^. 

Abundances are particularly interesting parameters since they can be measured quite di¬ 
rectly observing the Cosmic Microwave Background radiation. Accordingly, one has indepen¬ 
dent effective abundances HO (7 = 1,... ,n) which result from observations (the curvature HO is 
eliminated by the constraint X] ~ !)• More generally, a family of models is represented by 
a number of parameters e“ (the constants determining the function f{TZ)) and the independent 
abundances 

HO = H0(e«) (4.3) 


define a map from the set of models into the space of cosmic pies. 

By letting the parameters e“ vary, the equations (4.3) describe a parametrised subset S C M”. 
Each point on S represents the predictions of a model in the family. In this setting one can 
translate physical issues in geometrical claims, e.g., whether an observed cosmic pie (H?) G M” 
belongs to the image S C M"'. 

The subset S is called the model surface and it captures all the predictions of the models 
in terms of observable quantities. Thus if the observed cosmic pie is not on S then the family 
cannot account for observations and must be abandoned. 

In the model above, one has one parameter e and two independent abundances (H^,H))) G 
M^. The corresponding line representing the model is shown in Fig. 2. Since observations are 
considered to support (H^,H))) ~ (0.75,0.25) the family is discarded. 
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More generally, one can consider further observable quantities as the cosmographic parameters 
Id . Id . Id 

in the target space. On the other hand, ideally, one can expand the function fijV} in series 
(for example in Laurent series) and consider the infinitely many coefficients as the family pa¬ 
rameters e". Although, this becomes mathematically elaborated (having an infinite-dimensional 
surface embedded into an infinite-dimensional space) one can regard this as a universal setting 
to represent theories and observations as geometric intrinsic claims even when one may simplify 
it projecting to finite-dimensional representations in some particular context. For example, in 
this context is relatively easy and clear to argue if relations among observables hold in general 
or are valid within a specific family of theories when they emerge on the surface. 


5 Conclusions and perspectives 

In the literature. Palatini /(7^)-theories have not been studied much because of a number of 
shortcomings they are believed to have; see, for example, [8] and [25] where it was noticed that 
Palatini variation is not generally equivalent to metric variation. Since we totally agree on this 
point, we did not consider metric variations which are different theories. For us Palatini theories 
are motivated by the EPS formalism and they are, in our opinion, more natural than the metric 
ones. 

In the second place, one can perform a conformal on-shell transformation to show dynamical 
equivalence with purely metric Brans-Dicke theories (with a suitable potential); see [12, 15]. 
Since Brans-Dicke models (without a potential) are ruled out by solar system observations 
(e.g., precession of Mercury’s perihelia) it has been claimed that then any Palatini /(7^)-theory 
should be ruled out as well. This is false due to two issues. First, one should consider the role 
of the potential. Second, and more importantly, although field equations of the two theories 
are equivalent that does not imply a complete physical equivalence. In particular, in the case of 
precession of perihelia, the test on Brans-Dicke models is performed by assuming that planets 
go along timelike geodesics of g while in Palatini /(7^)-theories they go along geodesics of g. 
One can show by a detailed analysis that there is a huge family of Palatini /(7^)-theories which 
in fact pass the test of perihelia; see [19]. 

Moreover, the conformal transformation needed to recast Palatini-/(7^) theories into Brans- 
Dicke form is quite peculiar. One starts with a Palatini theory in which conformal transfor¬ 
mations act on the metric (which do not touch the connection which is independent) and they 
are symmetries. Then one solves the field equation related to the connection learning that the 
connection is in fact the Levi-Civita connection of a suitably defined conformal metric g. One 
is left with the field equation depending on the metric (and the conformal factor) alone since 
the connection has been eliminated. At this point by performing the inverse conformal trans¬ 
formation on the metric (which now affects the connection as well) the theory is recast into 
a Brans-Dicke form. There should be no surprise that the two theories are inequivalent since 
the latter transformation is not a symmetry. 

Another, issue is about star models. It was argued that assuming polytropic equation of state 
(EoS) one cannot match the inner solution with the outer solution for g. This is in fact true due 
to the fact that the conformal factor is continuous but not on the surface of the star. If one 
matches the conformal metric g (how is natural to do in view of the interpretation above) this 
is possible and the discontinuity appears only in the original metric g. Moreover, it has been 
noted that here the issue is due to the assumption about the EoS which is not a fundamental 
principle. Allowing a slight deformation of EoS one can in fact match the metrics at the surface 
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Figure 3. Approximation of a curvature density by cosmological constant and dust (dash line is the 
approximated density). 

of the star. The same result can be obtained by changing the external solution from a vacuum 
solution to an almost-vacuum one as it is reasonable and physically sound; see [3, 20, 21]. 

We showed how to classify extended cosmological models with respect to their predictions 
about the cosmic pie observations and to their predictions of the evolution of the scale factor. 

We showed that, the species revealed in the cosmic pie are essentially a choice in cosmology, or, 
in other words, an input from other physics (structure formation, baryogenesis, ...). However, 
once this input is given, extended gravity models produce quantitative predictions which allow 
in principle to falsify or corroborate models. 

The choice of effective species is not forced by mathematics essentially because of the appro¬ 
ximation method used. We chose to expand p in negative powers of a around the value a = 1. 
Even if we had p = 3a~^ and decide we want to see only dust and cosmological factor we would 
find 

2 

P ~ 1 + Vg, 

which approximates the actual effective density around today a = 1. 

One clearly sees that if one has a direct observation of a{t) in the past and assumes the 
approximated value then one can predict at what time in the past the deviation is first observed. 

Let us stress that, although it may not be computationally simple, the predictions of the 
models for the cosmic pie are extracted before solving the model. 

Future investigations are to be devoted to models in which the conformal factor stays positive 
everywhere, so that signature changes are prevented. 

Another research perspective will be investigating how one can encode physical properties 
(e.g., late time acceleration) as properties of the model surface and viceversa reading physical 
properties of the theories from the geometry of the surface. Currently, the main difficulty to be 
overcome is to control with enough generality the inversion of the master equation which plays 
a prominent role in expressing the embedding map (4.2). On the other hand, the model surface 
may be a good tool exactly because it shows that the properties are algebraically encoded in the 
model precisely through the master equation of the theory. 

A Conservation laws 

In a /(7^)-theory, the effective stress tensor is defined by equation (2.2). Its conservation 
with respect to the conformal metric g reads as V = 0, where indices have been raised by g. 
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The relations between controvariant stress tensors is then 

^ j,iiu _ 

Let us stress that indices of are here raised by g, while indices of are raised by g. 

Let us introduce the tensor 

AX = - {«}?;. = - 2X'5;)) V. Ill»), 

which accounts for the difference between the connections of the metrics g and g. Hereafter the 
covariant derivatives with respect to g will be denoted by while the ones with respect to g 

will be denoted by V^. The symbol means a covariant derivative when it is independent of 
any connection (in this case being applied to the scalar In Lp). 

One can easily show that 

In if 

and 

Then one has directly that 

= ln^(T^^ - f{n)g^^^- lTg^^^)gf,e - <7"“V„/(7^) - 

= + g"“(/"(7^)7^V,7^ - f{n)V^n - \V^T). 

However, by taking the covariant derivative of the master equation (2.1) one has 

/"(7^)7^V„7^ + /'(7^)V„7^ - 2/'(7^)Va7^ = ^V^T, 
so that one simply has 

so that is conserved with respect to g iff is conserved with respect to g. 

B Model solution 

Solving the FRW equations may be difficult though in a sense trivial since, in the end, one has 
to solve the Weierstrass equations with 

$(d) = -k + ^p(a(d))d^ 

i.e., computing the integral 



Maple can make it (with some effort) numerically (which means one can obtain a plot of 
a = a{t)). 
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Figure 4. Weierstrass function 4)(a). (Dash line is standard GR.) 


An easier way is to change variable back to {t, a), solve the equation and then compute the 
conformal factor, the new time and the conformal scale factor. This can be done by considering 
that 


da dd da dt 
di da dt di 




$(o). 


Since we know that di = y/(pdt and d{a) = a^Jdp one has for f{Tl) =Tl — 


$(a) := 


3 — €p{a) 
3 - epo 


dd . . 


-2 


$(p(a)). 


Once we have $( 0 ) we obtain the time coordinate as 

/ \ f°' da 

t{a) =to+ - 7 ^ 

Jao V^(a) 

and the conformal time as 

{t))dt = io+ j ^J^p{a)-^=. 

Jao V^(a) 

Finally the original density p is given as p{a) by EoS as well as the conformal density and 
pressure {p{a),p{a)). 

In our special case f{Tl) =TZ — setting e = | and po = 1 , one obtains (see Fig. 4) 

1 (2a3 - l)(5a3 - 

M (4a3 + 1)2 

with two allowed regions, one being a G [0, the other being a > 'pj^e first region cor¬ 

responds to an oscillating behaviour, hence to a closed universe. The second region corresponds 
to an expanding universe. Let us consider the second region. 

The original scale factor hence appears as a quite natural parameter to (more) easily integrate 
equations for d. 

In this way one is able to plot parametrically everything; see Fig. 5. Notice how the energy 
density of the visible matter p(t) stays bounded and the scale factor a{t) stays non-zero when 
standard GR has a singularity. 

This means that the solution can be extended to negative t and the big bang singularity is 
avoided (as some model of quantum gravity predicts; see [6]). 


i{a) = io+ 

Jto 
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Figure 5. a) Scale factor a as a function of t (dash line standard GR). b) Energy density p for visible 
matter as a function of t (dash line standard GR). 
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